Conformal Invariance and Shape-Dependent Conductance of Graphene Samples 
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For a sample of an arbitrary shape, the dependence of its conductance on the longitudinal and Hall 
conductivity is identical to that of a rectangle. We use analytic results for a conducting rectangle, 
combined with the semicircle model for transport coefficients, to study properties of the monolayer 
and bilayer graphene. A conductance plateau centered at the neutrality point, predicted for square 
geometry, is in agreement with recent experiments. For rectangular geometry, the conductance 
exhibits maxima at the densities of compressible quantum Hall states for wide samples, and minima 
for narrow samples. The positions and relative sizes of these features are different in the monolayer 
and bilayer cases, indicating that the conductance can be used as a tool for sample diagnostic. 



INTRODUCTION 

After the observation of the quantized Hall effect 
(QHE) in graphene [l|,lJ|, this material has quickly moved 
into the focus of research on quantum transport. Recent 
advances in patterning graphene into quantum dots [3|, 
ribbons [J, l5| and other shapes, and in selective gating 
of graphene sheets [6|, [Zl, |8|, has created a remarkably 
active field of graphene nanoelectronics. One of the chal- 
lenges in these experiments is finding reliable methods for 
sample characterization, allowing to distinguish graphene 
monolayer from graphene bilayer or multilayer systems. 
Quite often, the means for that are provided by Raman 
spectroscopy |9|. However, in a number of experiments 
it is more convenient to perform sample characterization 
using transport measurements [10| . 

The simplest transport characteristic of a graphene de- 
vice is its two-terminal conductance. This quantity ex- 
hibits plateau-like structure in the QHE regime, occur- 
ring at different electron densities in the monolayer and 
bilayer systems, which in principle makes it suitable for 
sample diagnostic. However, the observed conductance 
plateaus are often distorted, which is not surprising, be- 
cause the two-terminal conductance, unlike the resistiv- 
ity obtained from a four- probe measurement, in general 
depends on the sample aspect ratio and other geomet- 
ric characteristics. This dependence must be taken into 
account, in as much as possible, in interpreting the mea- 
surement results. 

The effect of sample shape on the conductance can be 
illustrated by the well known formula for a conducting 
square with the longitudinal and Hall conductivities Gxx 



and a^y and ideal contacts on opposite sides [11|, [1 



G 



L = W 



'xyi 



(1) 



which follows from a duality relation for 2d transport (see 
13, |l3[ and discussion below). The result ([T]) gives the 
macroscopic conductance in terms of microscopic trans- 
port coefficients. As we shall see, the dependence on 
axx and axy in Eq.|T]) is such that it can make the con- 



ductance Gl=w density-independent in the QHE regime 
near the graphene charge neutrality point (CNP), where 
both axx and (Jxy have strong density dependence. Thus 
the effect of sample geometry on conductance may be 
nonintuitive and should be accounted for carefully. 

Is it possible to extend the result ([1]), valid for a per- 
fect square, to other sample geometries? The next sim- 
plest shape to a square is a rectangle, for which the two- 
terminal conductance was found in Ref. Ill using confor- 
mal mappings of the Schwartz-Christoffel form. For all 
other shapes, luckily, the conduction problem can be re- 
duced, via a conformal mapping, to that of a rectangle 
(see Appendix). The aspect ratio of such an equivalent 
rectangle, which depends on the size of the contacts and 
on their separation, can serve as a parameter that classi- 
fies all conduction problems. 



In Ref.[ll[ a closed-form expression for the conduc- 
tance of a rectangle was obtained via an integral repre- 
sentation. However, as discussed below, the integrals of 
Ref. 11| are convergent very slowly, especially in the in- 
teresting limit of large Hall angles. Because of that, for 
our purpose it will be more convenient to employ the ap- 
proach of Rendell and Girvin lj| , which describes spatial 
distribution of the electric field and current in terms of 
a suitably chosen analytic function, allowing for direct 
numerical evaluation of the conductance. 

In this article, we use the method of Ref. 14] com- 
bined with the effective medium approach of Dykhne and 
Ruzin 15| that yields a semicircle relation between axx 



and axy We analyze conductance as a function of car- 
rier density, focusing on the features that distinguish be- 
tween transport in the monolayer and bilayer graphene. 
We conclude that the dependence of conductance on the 
sample shape, which may be quite strong, does not mask 
the difference between the monolayer and bilayer systems 
even in the absence of clear conductance plateaus. 



For the effective medium approach 15| to be appli- 
cable, the sample size must be large compared to the 
typical charge inhomogeneity length scale ^, otherwise 
strong mesoscopic sample-to-sample fluctuations are to 
be expected. In most of the paper we focus on the case 



of large samples, which can be described by a spatially 
uniform conductivity tensor obeying the semicircle rela- 
tion. We shall briefly discuss the situation in mesoscopic 
samples of size comparable to f at the end of the paper. 



DUALITY RELATION FOR CONDUCTANCE 

Here we shall focus on the rectangular geometry illus- 
trated in FigU] inset (later, in Appendix, it will be shown 
that for any conductor shape the problem can be mapped 
on that of an equivalent rectangle). To describe electric 
transport, we employ the bulk conduction approach, in 
which the sample bulk is characterized by the longitu- 
dinal and Hall conductivities Oxx-, <^xy The transport 
equation is j = ctE where ct is a 2 x 2 conductivity ten- 
sor, with the current and electric field obeying 



V-j = 0, VxE = 



(2) 



These equations must be solved with the boundary condi- 
tions jx = at X = 0, ly (current continuity) and E|| = 
at y = ±L/2 (ideal ohmic contacts). 

It is instructive to apply duality transformation [12lll3l| 
to this problem, rotating current j and electric field E by 
90° and interchanging them: j' = i?^/2E, E' ~ Rtt/2J- 
Upon such a transformation the transport equations ([2|). 
as well as the boundary conditions, preserve their form, 
whereby the conductivity tensor is replaced by a' = cr^^ 
and the dimensions of the rotated rectangle interchange: 
L' — W, W — L. Since resistance for the transformed 
problem R = V' /I' = 1/G" is identical to the con- 
ductance of the initial problem G = I/V, where V is 
source/drain voltage and / is net current, we obtain a 
duality relation 



G{L,W,a) = G-\L',W',a') 



(3) 



We note that, since a' = a^^ is the resistivity tensor, the 
quantity G(L', W , a') has dimension of resistivity, and so 
the right hand side of ^ has dimension of conductance. 
To simplify the relation ([3]), we take into account that G 
scales with ct, i.e. that G{L,W,ri(j) = riG{L,W,a), and 
that it is invariant upon sign reversal of a^y Writing 



<7xx/i<7. 



XX ' xy ) 



(^'xy = (^xy/{(jlx+'^xv), ^ud USing 



the scaling property of G, we obtain a relation 



G{L,W,a)^{ai,+aiy)/G{W,L,a) 



(4) 



which connects the rectangles L x W and W x L having 
the same bulk transport coefficients. Setting L = W, we 
obtain the conductance of a square, Eq.([T]). 



THE SEMICIRCLE MODEL 

The effective medium approach [l5'| provides a con- 
venient framework for understanding the density depen- 
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FIG. 1: Longitudinal and Hall conductivity for (a) graphene 
monolayer and (b) graphene bilayer, obtained from the semi- 
circle model, Eqs.®,!!?!,®, for two values of the Landau 
level width parameter A = 1.7 (solid lines), A = 0.5 (dashed 
lines). Inset: Schematic of a conducting sample of dimensions 
L and W, of a rectangular shape, with source and drain at 
opposite sides. 



dence of axx and axy in QHE systems. It predicts a semi- 
circle relation between Gxx and Gxy, derived from a two- 
phase model, in which the system at the QHE transition 
is treated as a mixture of incompressible puddles with lo- 
cal Hall conductivities a' , a'', given by the quantized 
values at the neighboring QHE plateaus. 

The semicircle relation is a statement about proper- 
ties of the macroscopic conductivity tensor on the length 
scales much greater than the size of individual puddles. 
Although the validity of this relation, strictly speaking, 
is limited to the regime dominated by large-scale fluctu- 
ations of electron density, it is empirically known to pro- 
vide a good description of the integer QHE plateau tran- 
sitions observed in various semiconducting systems [iq . 
Recently, the semicircle relation was employed to de- 
scribe transport coefficients in graphene ITJ, [iSl, ll9| . 

We stress that, while the semicircle model is realistic. 



and also quite convenient to use, its specifics are almost 
certainly not essential for our conclusions. We believe 
that slight departure in behavior of transport coefficients 
from the semicircle model will have little effect on the 
properties of conductance. 

Prior to turning to the semicircle model, we recall 
that the conductivity a^y in graphene monolayer ex- 
hibits steps of size 4e^//i between adjacent integer quan- 
tum Hall plateaus, where the factor of four describes 
combined spin and valley degeneracy of Landau levels. 
The incompressible densities corresponding to the QHE 
plateaus in graphene monolayer are [ij, |j| 



z/„ = 4(n+l/2)|B|/$o, n = 0,±l,±2. 



(5) 



where |i?|/$o is electron density for a single Landau level. 
In graphene bilayer QHE, due to accidental degeneracy 
of the Landau level positioned at the neutrality point h' = 
0, there is an 8e^/h Hall conductance step between the 
plateaus with axy = ±4e^//i, whereas other conductance 
steps are of normal 4e^//i size. Accordingly, in the bilayer 
the incompressible QHE densities are j20j 



4n|B|/$o, 



±1,±2. 



(6) 



In both cases ([5]) and ([6]) the density values i^n are ar- 
ranged symmetrically around the neutrality point. Hall 
conductivity on the plateaus takes the values (Jxy\n = 
Vnc'^/h, where Vn are densities ^ and ([6]). 

In the semicircle model 15], the contributions of each 
Landau level to the longitudinal and Hall conductivities 
5n'^xx{v), SnCxyi^^) are related by the semicircle law: 



^nCTxx + {SnCTxy - 4y,n) ('^riCT^ry 



.(0) A- 



= 0, (7) 



where i7iy,n and a^^ „, are the quantized Hall conductiv- 
ities on adjacent plateaus. Here n and n' are neighboring 
integers in the sequence ... — 2,— 1,0,1,2... for the mono- 
layer, and ... — 2, —1, 1, 2... for the bilayer: n' = n + 1 
except the double-degenerate v ~Q Landau level for the 
bilayer, in which case n = —1, n' ~ 1. 

The longitudinal conductivity (J,iCr^a;(i^) exhibits a peak 
centered at the Landau level. We model it by a gaussian 



OnOxx{v) = -C^e ^ 2V"T„;^ 



(8) 



where the parameter A describes broadening of the Lan- 
dau level (large values of A correspond to a narrow Lan- 
dau level). 

In the semicircular model, the peak value of ^n^xx 
must equal to ^{crly n' ~ ^xy,n)- This is ensured by the 
prefactor C„ value in Eq.® chosen to coincide with the 
nth Landau level degeneracy. For graphene monolayer 
we set C„ = 4 (spin and valley degeneracy) for all n, 
while for the bilayer C„ = 8 for n = — 1 (spin, valley and 
accidental degeneracy) and C„ = 4 for all other n's. 



The total conductivity tensor is given by the sum of 
the contributions of all Landau levels, 



(Jxx{l^) = '^Snaxxil'), CTxyW) = '^SnCTxyi^)- (9) 



For simplicity, here we choose the same value of the pa- 
rameter A for all Landau levels. The resulting conduc- 
tivity density dependence is illustrated in Fig[T^,b. 

We point out that an interesting prediction can be 
drawn, specific to the graphene zeroth Landau level 
{v — 0), by combining the semicircle model ([7]) with 
Eq.([T]). In a square sample with a negligible overlap be- 
tween Landau levels, the two-terminal conductance ([T]) 
would be completely density-independent across the ze- 
roth Landau level, because in this case \ a^x + ^xy would 

2 2 

equal 2^ for the monolayer and 4^ for the bilayer. This 
happens because the density dependence of transport co- 
efficients, the peak in axx and the step in axy, centered at 
I' = 0, cancel each other in the expression ^. For weakly 
overlapping Landau levels, the contributions of the levels 
adjacent to z^ = would lead to slight deviations from a 
fiat plateau. 

The conductance measured in graphene indeed often 
exhibits a plateau across the entire v = Q region. Two ex- 
amples of such behavior in recent literature are Ref . [22] , 
Fig. Id, and Ref. 7], Fig. 3B. In both cases, the measured 
conductance is nearly flat in a wide density interval cen- 
tered at I' = 0, with a small peak in the middle. Interest- 
ingly, the sample geometry in both cases was quite differ- 
ent from a square. In Ref. [22] the AFM image indicated 
that the sample was approximately rectangular with con- 
tacts at opposite sides, similar to the schematic in Fig[T] 
inset. However, its width was quite large, W ~ 5L, in 
which case a fairly large peak at z/ = is to be expected 
(see below). To understand the discrepancy, it would 
be useful to known how spatially uniform the conduc- 
tion was, in particular near contacts (poor contact along 
part of the sample edge could reduce the effective sam- 
ple width, bringing W closer to L). As for the device 
in Ref. 7], its contact geometry was quite far from rect- 
angular, yet giving rise to a fairly good plateau around 
i/ = 0. 

Such a behavior can be understood in the context of 
universality of the conductance problem that results from 
its conformal invariance. As we discuss in Appendix, for 
an arbitrary geometry of a conductor with two contacts 
of arbitrary sizes the conductance is identical to that of 
an equivalent rectangle. In other words, all possible con- 
duction problems are classified by the values of a single 
parameter, the aspect ratio L/W of an equivalent rectan- 
gle. Then, as long as L/W « 1, the conductance would 
behave in the same way as for a square shape, even if the 
actual geometry is very different from a square. 
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FIG. 2: Two-terminal conductance (I15|l of a rectangular 
graphene sample: (a) monolayer, (b) bilayer, for aspect ra- 
tios L/W — 0.25, 0.5, 1, 2, 4 (top to bottom) with the Landau 
level width parameter A = 1.7 (corresponding to solid lines 
in Fig[l|. Arrows mark the incompressible densities ([5|l. (pjl. 
Note the plateau at ;/ = for the square case, L — W (red 
curve), which is in agreement with the behavior predicted by 
Eq.©. 



CONFORMAL MAPPING APPROACH 

Now we proceed to describe the distribution of cur- 
rent and electric field in a rectangular sample with an 
arbitrary aspect ratio L/W and spatially uniform con- 
ductivities axx, <^xy (see Figlljinset). This problem, with 
zero boundary condition for normal current component 
at the sample boundary and for tangential field compo- 
nent at the sample/contact interface, has been treated in 



Refs. [111 . |14| using conformal mapping technique. 



In Ref. [11| the Schwartz-Christoffcl mapping was em- 
ployed, leading to the following integral representation 
for the conductance: 



G 



■ J(l,-1) 

/(iA,i) 



(10) 



where 



/(m,u) 



di 



|(C-i)(iA + OI*+I(C + i)(iA-OI'-' 

(11) 

Here 8± = 1/2zL9/tt, and 9 = tan ^ {ayx / (Txx) is the Hall 
angle. The parameter k above is the elliptic modulus, 
< k < 1, related to the sample aspect ratio via 



L 
W 



K{k) 
2K{k') 



k' = vT^^, 



(12) 



where K{k) is the complete elliptic integral of the first 
kind. 

In principle, Eas. pU)) . pT|) . p^ give a complete solu- 
tion of our problem. However, the integrals in Eq. ljlOp 
contain power-law singularities with the exponents 5± = 
1/2 ± 9/tt that can approach unity for 9 w ±7r/2. This 
makes the integrals in Eq. ()11|) difficult to evaluate numer- 
ically, especially for large Hall angles, when the singular- 
ities are the strongest and the integrals are converging 
fairly slowly. 

Instead, we use a different approach, developed in 
Ref. [iJ] , which provides an expression of the electric field 
and current in the sample in terms of the exponentials 
of infinite but rapidly convergent sums. This method is 
more convenient for our purposes, because it allows to 
choose the integration contour for numerical evaluation 
of the conductance so that it bypasses the singularities. 

In Ref. [14| , the electric field components Ex and Ey at 
a point z — x + iy are obtained as the real and imaginary 
parts of a suitable analytic function. The latter function 
is found to be of the form 



Ey + iEx = -e 



fi^) 



where 



f{z) = ^9- Y. 



n>0(odd) 



49 smh{mTiz/W) 
nil cosh(n7rL/2W^) 



(13) 



(14) 



Current distribution can be obtained by combining (|13p 
with the relation j^. -|- ijy = {oxx + io-yx){E,^ + iEy). 

Once the current distribution is found, it can be used 
to obtain the two-terminal conductance 



G = I/V, 



(15) 



where / is the total current and V is the source-drain 
bias voltage. To evaluate the net current / by integrating 
current density, one has to choose a cross-section through 
the sample that does not pass through its corners, where 
the function f{z) has singularities. It is particularly con- 
venient to perform this integration along a straight line 
that cuts through the middle of the rectangle at y = 
(along the x-axis in Fig[l] inset): 



/ = 



w 



jy{x,0)dx. 



(16) 



Since jy = (T^xEy 



UyxEx, the current is given by 



w 



c{Ey{x,0) +tSin9Ex{x,0))dx. 



(17) 



For similar reasons, we calculate the voltage drop be- 
tween the upper and lower contacts (y = ±Z//2) as an 
integral of the electric field along a straight line con- 
necting the points (VF/2, -i/2) and {W/2, L/2) of the 
contacts: 



V 



L/2 



L/2 



Ey{W/2,y)dy. 



(18) 



Evaluating the integrals (fT7| and (|18p numerically, we 
obtain the conductance (|15p as a function of transport 
coefficients, which defines its dependence on electron den- 
sity. The results for different aspect ratios L/W for the 
monolayer and the bilayer case are displayed in Fig[2]and 
FiglHl for two different values of the Landau level width 
parameters in ([8]), A — 1.7, 0.5. 

As evident from Figlll the conductance behavior de- 
pends strongly on the sample aspect ratio L/W. For a 
square sample, L = W, the conductance is a monotonic 
function of density at positive and negative ;/, rising in 
steps associated with QHE plateaus. In this case, the 
behavior of G is found to be completely flat near 1^ = 
in agreement with the above discussion based on Eq.([l]). 

For L 7^ W, the conductance turns into a nonmono- 
tonic function of density, with the QHE plateaus being 
less pronounced than for L = W^. For L < W the conduc- 
tance has minima near QHE plateau centers, Eq.((5|) for 
the monolayer, and Eq.® for the bilayer case (marked 
by arrows in Figs l2|3p . while for L > W^ the conductance 
has maxima at these densities. Overall, the conductance 
behavior ioi L < W is "inverted" compared to that for 
L > W, as expected from the duality relation ([4]) which 
implies 



G 



L<W 



> Gl—\y > G 



L>W 



for all I'. Recently, the conductance as a function of 
carrier density was studied for samples of several aspect 
ratios [2ll | , with the results being in qualitative agreement 
with ours. 

With these results at hand we are now in position to 
ask which differences between the monolayer and bilayer 
systems are robust with respect to Landau levels' broad- 
ening and variation in sample geometry. First, we note 
that the centers of the QHE plateaus (arrows in Figs l2|3|) 
remain close to the positions of the minima in G for 
L < W, and its maxima for L > W. The corresponding 
incompressible values of the density, Eqs.®,®, which 
are marked by arrows in Figsl2l[3l are equally spaced in 
the monolayer case, but are not equally spaced in the bi- 
layer case due to the eight-fold Landau level degeneracy 
at 1/ = 0. It can be seen by comparing Fig[2]and FigO 
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FIG. 3; Same as in Fig[2]for broader Landau levels, described 
by ([SI with A — 0.5 (dashed lines in Fig[l}; the sample aspect 
ratios are L/W = 0.25, 0.5, 1, 2, 4 (top to bottom). Note that 
the qualitative features, such as the positions of the conduc- 
tance minima at the QHE plateau centers for L < W (max- 
ima for L > W), as well as the conductance values at these 
densities, are similar to those seen in Fig[2] despite increased 
Landau level broadening. Note also the relative size of the 
1/ — peak in the monolayer and bilayer cases, compared to 
the size of neighboring peaks at other compressible densities, 
which is also insensitive to the Landau level broadening. 



that this difference between the monolayer and bilayer 
systems is not masked by Landau levels' broadening or 
by the aspect ratio variation. 

Next, the conductance values at the minima for L < 
W and at the maxima for L > W remain close to the 
associated QHE values G « VnC^ /h, which are different in 
the monolayer and bilayer cases. This difference is clear 
in Figs l2l3l even for large aspect ratios and broadened 
Landau levels. In practice, however, these values may 
change as a result of added contact resistance. If this is 
the case, the relative positions of the plateaus in density, 
inferred from the arrangement of maxima and minima of 
G, are more robust than the conductance values at these 



densities. 

Finally, we note the difference in the size of the peak 
at zy = as compared to the sizes of neighboring peaks 
at other compressible densities, found for the monolayer 
and bilayer at L < W. (The same is true for the dip at 
ly = ior L > W.) In the monolayer case, the size of 
the 1/ = peak/dip is comparable to the sizes of neigh- 
boring peaks/dips, whereas in the bilayer case the u = 
peak/dip is almost two times larger than its neighbors. 
This difference reflects the higher density of states in the 
1/ = Landau level of a bilayer. 

It is also interesting to note that the heights of all peaks 
and dips in G are completely independent of the Landau 
level broadening, as can be seen by comparing the curves 
in Figl2] and FiglS] This behavior is specific to the semi- 
circle model, in which the peak values of axx and corre- 
sponding values of axy are universal. Indeed, as evident 
from FiglH peaks of axx line up with steps in axy, and the 
values of axy in the middle of each step are independent 
of the Landau level broadening. As a result, although 
broadened Landau levels change the overall behavior of 
the conductance G, the peak values remain intact. 




FIG. 4: Examples of possible p-type and n-type puddle config- 
urations near graphene charge neutrality point for mesoscopic 
samples with at least one of the dimensions comparable to the 
typical puddle size 5. In (a) and (b) two possible configura- 
tions are shown for a square sample with L, W ^ £,. The 
corresponding conductance values, G = ie^/h for (a) and 
G — e^/h for (b), are different from the value G* = 2e^ /h 
predicted by the effective medium (semicircle) model. Pud- 
dle arrangement in a short and wide sample {L '^ ^ <^ W) 
and in a narrow and long sample {W ^ ^ <€. L) is illustrated 
in (c) and (d). 



CONDUCTANCE FLUCTUATIONS IN SMALL 
SAMPLES 

So far we have described the system by a spatially uni- 
form conductivity, assuming this effective medium model 
to be valid at all filling factors. This assumption, ade- 
quate for large samples, breaks down for smaller samples 
when the sample dimensions L or W become compara- 
ble to the typical charge inhomogeneity length scale ^. 
Conductance of such small samples at QHE transitions 
can exhibit strong sample-to-sample fluctuations and sig- 
nificant deviations from the prediction of the semicircle 
model. 

The effect of charge inhomogeneity is expected to be 
especially strong near the charge-neutrality point (CNF), 
where the density of free carriers is low and screening of 
the disorder potential is poor, for both the monolayer and 
bilayer graphene. For simplicity, below we shall focus on 
the monolayer graphene in the vicinity of the CNP, i/ « 0. 
To qualitatively understand the effect of charge inhomo- 



geneities, we employ the two-phase model [15[ treating 
the sample as a mixture of incompressible puddles of 
types n and p of typical size ^. The filling factors in the 
puddles are ly — ±2, corresponding to the QHE plateaus 
adjacent to CNP. 

To illustrate the effect of strong fluctuations near the 
CNP, we first inspect the case of a small square sample, 
with both dimensions L and W comparable to £,. Suppose 
that the carriers in such a sample form just two puddles, 
of type p and n (j/ = ±2). There are two topologically 
distinct puddle configurations giving the net conductance 
different from the value G* = 2e^/h expected from the 



semicircle relation: (i) the p-n boundary connecting the 
source and the drain (see Figd^) (ii) the p-n boundary 
connecting the opposite free edges of the sample (see 

FigHb). 

In the case (i) the p and n regions are connected in 
parallel, and thus the net conductance is the sum of the 
conductances Go = 2e^//i, giving G — 2G* — ie'^/h. 
In the case (ii) the p and n regions are connected in se- 
ries, with the net conductance equal G = 1/2G* = e^//i 
(see Ref. [23[). Thus for mesoscopic square samples with 
L, W '^ (, the conductance strongly depends on the pud- 
dle configuration and can exhibit either a peak or a dip at 
1^ = 0, whereas the semicircle model predicts a plateau. 

Now let us consider a sample which is short but wide: 
L < S, -^ W (see Fig|3t). Such a sample consists of 
about N ^ W/^ alternating p and n regions of width 
~ ^. The total conductance is given by the sum of the 
conductances of the individual regions, the conductance 
of each region being Go — 2e'^ /h. This yields the net 
conductance of G - NGq = {2W/0^- This is about 
i/^ smaller than the value G* = {2W/L)^ predicted by 
the semicircle relation for the effective medium model. 

Similarly, for the case of a narrow but long sample, 
W ^ ^ <^ L (see FigHJi) the conductance is given by the 
series conductance of A^ ~ W/^ puddles: G ~ (2^/L)^. 
This is £,/W times larger than than the prediction of the 
effective medium model. These estimates show that near 
CNP the semicircle model overestimates the conductance 
of very short samples, and underestimates the conduc- 
tance of very narrow samples, when at least one of the 
sample dimensions is comparable to the puddle size. 



CONCLUSIONS 



(a) 



(b) 



(c) 



The above results for the two-terminal conductance 
of rectangular samples, and in particular, the conduc- 
tance dependence on the sample aspect ratio, can serve 
as a benchmark for understanding properties of graphene 
samples. We found that G exhibits peaks (dips) at the 
compressible densities for L < W {L > W), which com- 
pletely disappear at L — W. We could identify several 
specific features that may help to distinguish between 
transport in graphene monolayer and bilayer. Those in- 
clude positions of the incompressible densities, inferred 
from minima (maxima) in G a.t L < W {L > W), the 
values of G at these densities, and the relative size of the 
central peak (dip) in G as compared to the neighboring 
peaks. These features, which are shown to be insensitive 
to the sample aspect ratio and to Landau levels' broaden- 
ing within our model, can be used for sample diagnostic 
in transport measurements. 
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APPENDIX: RECTANGLE AS THE MOTHER OF 

ALL SHAPES; CONFORMAL INVARIANCE AND 

UNIVERSALITY 

There is a profound relation between the two-terminal 
conductance and conformal invariance of the 2d trans- 
port problem. It arises because for 2d conductors of arbi- 
trary shapes with spatially uniform axx and axy the con- 
ductance is invariant under conformal transformations, 
and because all single-connected domains in the plane 
can be conformally mapped to one another. As a result, 
a conducting domain of any shape has the same con- 
ductance as a particular domain of some simple shape 
with an appropriate arrangement of contacts. The sim- 
ple shape can be chosen in a number of ways, in par- 
ticular it can be chosen to be a rectangle. We show in 
this section that for a conductor of any shape with any 
configuration of two contacts the conductance is equal to 
that of a rectangle with some aspect ratio L/W . 

Because the correspondence between conductors of ar- 
bitrary shapes and equivalent rectangles is purely geo- 
metric (it is defined by a conformal mapping) , the aspect 
ratio L/W of an equivalent rectangle depends on the sam- 
ple shape but does not depend on the values of transport 
coefficients Gxx and (Jxy As a result, one can use the 
rectangle problem with a fixed aspect ratio to describe 
conductance as a function of the carrier density, via the 
axx and axy density dependence. 

To formulate the constraints due to confornral invari- 
ance, we recall that conformal mappings in 2d are realized 
as analytic functions of the complex variable z = x -\- iy. 




FIG. 5: Conformal invariance illustrated by several conduct- 
ing domains with contacts. If two domains can be conformally 
mapped on each other so that the contact regions are mapped 
on the corresponding contact regions, their conductances are 
the same. The Riemann mapping theorem guarantees exis- 
tence of a mapping between an arbitrary domain (a) and a 
unit disk (b) such that three points on the boundary of (a), 
marked 1, 2 and 3, are mapped on any three points (1,2 and 
3) on the circle (b). The position of the fourth point, which is 
not specified, defines a one-parameter family of possible con- 
duction problems. All such problems can be parameterized 
by rectangles (c) with difi'erent aspect ratios. 



Thus we consider mappings between two complex planes 
z and w defined by analytic functions w ~ f{z), which 
map the sample domain (hereafter denoted V) in the 
plane z onto a domain V in the plane w. 

The relations ^ as well as the boundary conditions 
-E|| = on the contacts and j± = on the open 
boundary are invariant under such conformal mappings. 
The easiest way to verify this is to note that the cur- 
rent continuity condition V • j = can be solved by 
j = V X {zip{x,y)), where z is the unit vector normal to 
the plane. Then the relation V x E = combined with 
ix + ijy = io'xx + io'yx){Ex+iEy) means that the function 
if is harmonic, i.e. it satisfies the Laplace's equation 



{dl 



dl)^{x,y) = 0. 



(19) 



Because a harmonic function remains harmonic under a 
conformal mapping, and because the angles between the 
gradient Vi^(z) and the boundary of T) are the same as 
the angles between W(p{w) and the boundary of V at 
corresponding points, Eqs.© as well the relations _B[| = 
and j±_ = are indeed conformally invariant. 

This simple mathematical fact can serve as a platform 
for the following interesting observation. Suppose an an- 
alytic function w = f{z) maps the conducting domain T> 
in the plane z to a domain T)' in the plane w so that a 
pair of contacts to V is mapped on a pair of contacts to 
V . More precisely, let four points 1, 2, 3, 4 on the bound- 
ary of V be mapped onto four points l',2',3',4' on the 
boundary of V . Let the arcs 1 — 2, 3 — 4 and 1' — 2', 
3' — 4' be ohmic contacts for the problems in the z and w 
planes, respectively (see Fig (5]). Then the two-terminal 
conductance of 2? is exactly the same as that of 2?', i.e the 
conductance is invariant under all conformal mappings 
that map to one another the corresponding conducting 
domains and contacts. 



t y 




the same: 



FIG. 6: Complex halfplane and conduction problem in it, 
which can be mapped on that in a rectangle using the 
Schwartz-Christoffel mapping, Eq.((23}. The aspect ratio 
L/W of the equivalent rectangle depends on the relative size 
of the contacts, shown in blue, and the distance between them. 
The end points of the contacts are ^1,2 = Tli ^3,4 ~ ±l/fc. 



On the other hand, as is well known from the theory of 
complex variables (the Riemann mapping theorem), any 
two single-connected domains V and V can be mapped 
onto each other. This mapping can be fixed so that any 
three points on the boundary of V are mapped to any 
three points on the boundary of V (under these condi- 
tions the mapping is unique). This allows to reduce the 
conductance problem of an arbitrary domain V to that 
of some simple domain, e.g. a circular disk \w\ < 1. 

Thus the only reason the conductors of different shapes 
do not all have the same conductance is the additional 
freedom in choosing the contacts, defined by the points 
1,2,3,4. Furthermore, because any conduction prob- 
lem is equivalent, via a conformal mapping, to a disk 
with contacts defined by four points on the boundary 
zi, ..., Z4 — e*^S ...,e*^*, and because three of those points 
can be fixed by the Riemann mapping theorem, the only 
freedom left is in the position of the fourth point. There- 
fore, all conductance problems form a one-parameter 
family. For the points zi,...,Z4 this parameter can be 
expressed, e.g., as the so-called cross ratio 



A1234 — 



(zi - Z4){Z3 -Z2) 
(zi - Z2)(Z3 - Z4) 



(20) 



which takes real values for any four points that lie on a 
circle. Applied to the problem of a rectangle, this pro- 
cedure would give a one-to-one relation between the pa- 
rameters L/W and A1234, proving that indeed any con- 
ductance problem is isomorphic to that of a rectangle. 

To obtain such a relation between L/W and A1234, let 
us consider a complex halfplane Im ^ > with four points 
on the real axis 



a,2 = Ti, e3,4 = ±i/fc, 



(21) 



with the parameter < A: < 1 (see Figl6]). Under a 
fractional-linear mapping that maps the unit disk to the 
halfplane, the points zi, ...,Z4 can be mapped onto cor- 
responding points (|2ip if and only if the cross ratios are 



^1234 



(1 - fc2)/4fc 



(22) 



(the cross ratio is an invariant of fractional linear 
mappings). On the other hand, using the Schwartz- 
Christoffel formula, the halfplane Im^ > can be 
mapped on a rectangular domain by the function 



Jo 



de 



v/(i-e'2)(i-fc2e'2)' 



(23) 



where F{£^, k) is the elliptic integral of the first kind. This 
function maps ^^2 — Tl to zi_2 — TK{k), where 



K{k) 



de 



^{i-e){i-k^e) 



(24) 



is the complete elliptic integral of the first kind. The 
points a, 4 = ±l/fc are mapped to Z3,4 = zLK{k)+iK{k'), 
where k' = Vl — k^. Thus the sides of the rectangle have 
lengths W = 2K{k) and L = K{k'), which gives the re- 
lation between L/W and k of the form p^ . Combin- 
ing this with (|22|) . we can relate L/W to the cross ratio 

Ai234- 

In the halfplane of Figl6]the distribution of the electric 
field and current can be found by noting that the function 
E{£^) = Ey + iEx is analytic at Im^ > and at the 
boundary its argument takes fixed values between the 
points a,. ..,^1 argS_i<^<i = 0, arg£'i<5<i/fe ^ f + 6', 
argS_i/fc<^<_i = -|-f e, arg£:|5|>i/fe = -tt. These 
requirements are sufficient to reconstruct the function: 



E{0 



(i-e)*+(i + o'-(i-fcC)'-(i + fcO'+ 



(25) 



where 5± — ^if with 9 = tan^^ Oyx/o^^ the Hall angle, 
as in (jlip , and the unknown prefactor A depends on the 
source-drain voltage. Current distribution is then found 
as J(^) = jy + ijx — [oxx + i<^xy)E{(/). The total current 
and voltage can be expressed by integrating J(^) over the 
contour C/, and E{£/) over the contour Cy (see Figl6]). 
The ratio of these integrals gives the result (fTO|) . 

To summarize our discussion, because the conductor 
of an arbitrary shape can be conformally mapped on a 
rectangle which has the same conductance, the rectangle 
problem is "universal." For a sample of any shape with a 
pair of contacts of arbitrary size and form, an equivalent 
rectangle can be found such that it has the same conduc- 
tance. Crucially, the aspect ratio L/W of this rectangle is 
independent of the values of transport coefficients, which 
means that its conductance will have the same density 
dependence as that of the physical sample. 

We note that, although the results of this section ap- 
ply to conductors of completely arbitrary shapes, and in 
that sense they are far reaching, there are several limita- 
tions. First, we have assumed that the electron system 



is spatially uniform and homogeneous, i.e. the trans- 
port coefficients Uxx and a^y are position- independent. 
Our second assumption was that transport is fully de- 
scribed by the 2d current-field relation j = ctE, where a 
is the conductivity tensor. In particular, the above model 
does not allow for the edge current states, which could, 
in principle, alter the shape dependence of the conduc- 
tance. Still, since the bulk transport model agrees with 
the edge transport model in the limit of a large Hall an- 
gle, it can probably provide a good guidance even when 
the quantum Hall effect is fully developed. 
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